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Spinor Slow-Light and Dirac particles with variable mass 
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We consider the interaction of two weak probe fields of light with an atomic ensemble coherently 
driven by two pairs of standing wave laser fields in a tripod-type linkage scheme. The system is 
shown to exhibit a Dirac-like spectrum for light-matter quasi-particles with multiple dark-states, 
termed spinor slow- light polaritons (SSP). They posses an "effective speed of light" given by the 
group-velocity of slow-light, and can be made massive by inducing a small two-photon detuning. 
Control of the two-photon detuning can be used to locally vary the mass including a sign fiip. This 
allows e.g. the implementation of the random-mass Dirac model for which localized zero-energy 
(mid-gap) states exist with unsual long-range correlations. 
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Due to the recent experimental and theoretical 
progress on various system there has been a growing in- 
terest in interacting systems with effective single-particle 
Dirac dynamics. A prime example is Graphene [1] which 
shows unusual low-energy properties, described by Dirac 
quasi-particles. The effective "speed of light" of these 
quasi-particles is given by the electron Fermi velocity, 
which is about 300 times smaller than the vacuum speed 
of light. We here show that an effective Dirac-like dy- 
namics also emerges for two-component, i.e. spinor-like, 
slow-light polaritons [2] in ID. They possess a control- 
lable effective speed, corresponding to the much smaller 
group velocity of light [3] in media exhibiting electromag- 
netically induced transparency (EIT)[2]. This may allow 
experimental studies of a number of interesting effects of 
relativistic dynamics at low energies and small velocities. 
Additionally the effective mass, determined by laser de- 
tuning, is locally and dynamically adjustable including 
the sign, which provides access to interesting phenomena 
such as the unusual localization in a Dirac model where 
the mass is a random function of space [4, 5]. The lat- 
ter is inaccessible for truly massive particles such as cold 
atoms. 

Slow-light polaritons are formed in the Raman inter- 
action of a weak probe field with a coherently driven en- 
semble of atoms with a A- type linkage pattern (Fig. la). 
They build the basis of ultra slow light [3] , light storage 
[6, 7] and, upon using two counter-propagating control 
fields (Fig. lb, c), stationary light [8-10]. They exhibit ex- 
traordinary long lifetimes ranging from several hundreds 
of microseconds [6] up to hundreds of millisecond [11, 12] 
and offer a wide range of tuning parameters. 
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FIG. 1. (Color online) (a) 3-level A-type linkage pattern for 
the creation of slow-light polaritons. The strong control field 
Q produces electromagnetically induced transparency (EIT) 
for the weak probe field £. (b) Linkage pattern with two 
counter-propagating control fields Q± to create a stationary 
pattern of counter-propagating probe light £±. (c) Typical 
experimental setup for the creation of single-component sta- 
tionary light by use of two counter-propagating, mutually or- 
thogonal control fields Q±. (d) Tripod-type linkage pattern 
for the creation of spinor slow-light polaritons (SSP). 



In the following we argue that by adding additional 
states to the above linkage pattern (Fig. Id), we can cre- 
ate a spinor-like object consisting of two adiabatic eigen- 
solutions immune to spontaneous decay. The suggested 
tripod linkage interaction between atoms and lights is a 
minimal realization of spinor slow-light polaritons (SSP). 
They obey an effective ID Dirac equation, thus creat- 
ing relativistic quasi-particles with an internal degree 
of freedom. They travel at an effective speed c* given 
by the slow-light group velocity, and posses a mass m* 
that is determined by a variable two-photon detuning. 
The possibility of a locally adjustable mass, which is ab- 
sent for any truly massive particle, allows to study a 
number of interesting phenomena: For example, if the 
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FIG. 2. (Color online) Dispersion relation of the double- 
tripod pattern (Fig. Id) (ignoring decay) showing the energy 
branches of the spinor slow-light polaritons (SSP) (red dotted 
line: 5 = 0, blue solid line; S = 0.075^0 ). 



mass of the Dirac particle is a randomly varying func- 
tion of space with a vanishing mean-value, there exists 
a mid-gap (zero-energy) state with unusual correlations. 
Random-mass Dirac Hamiltonians describe a number of 
effects in condensed-matter systems, ranging from disor- 
dered half- filled metals [13], random anti- ferromagnetic 
spin-1/2 chains [4, 14, 15], random transverse-Ising spin- 
1/2 chains [16, 17], spin Peierls and spin-ladder systems 
[5, 18, 19]. Many aspects of random-mass Dirac systems 
arc not understood, e.g. the effect of local impurities or 
boundaries [20] or when interactions are added. Spinor- 
polaritons may offer new experimental access to these 
issues. 

We here consider light coupled to an ensemble of atoms 
with a linkage pattern as seen in Fig. Id, consisting 
of three (meta-)stable states \g), \s), \h) and one ex- 
cited state |e) with a decay rate 7. There are two 
pairs of counter-propagating control fields with Rabi- 
frequencies ri^j-e^^'^^^ and fi^.-i-e^'*''''^ on the transitions 
|s) — |e) and \h) — |e), respectively, as indicated in Fig. 
Id. These create KIT for a pair of counter propagat- 
ing probe fields E±, coupling the transition \g) — |e) [2]. 
The use of pairs of counter-propagating control beams 
is essential to produce two-component slow- or station- 
ary light, a feature missing in the ordinary tripod-type 
linkage pattern [21-23]. This scheme effectively creates 
two parallel tripod-type linkage patterns sharing com- 
mon ground-states. We introduce amplitudes £±{z,t) of 
the probe fields that vary slowly in space and time by 
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£±{z, t) eyi'p{—i{ujpt^ kpz)} + h.a. . Further- 
more we define continuous flip operators via d^i,{z^t) = 

-KnT.j(^av{z)^U^ ^i*l^ ^U = Ia^)m(''I bemg the flip 
operator for the j'-th atom and the sum being taken over 
AiV atoms in a volume AV{z). 

The full dynamics of the system is governed by the 
Maxwell-Bloch equations (MBE) in ID. In the weak field 
limit the ground-state is only slightly depleted allow- 
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ing to assume agg « 1 and a^^ « in lowest order of 
perturbation. We further expand the relevant coherence 
CTge = (Tge e^'^^ + crgee~^^'^, and neglect all higher spa- 
tial fc-components as well as couplings to other coher- 
ences with higher Fourier /c-modes. This approximation 
is know as the secular approximation [8, 24], and is well 
justified in hot atomic gases [9, 25]. We can thus write 
down the linearized Maxwell-Bloch equation (MBE) in 
form of a 6 X 6 matrix consisting of 2 x 2 sub-blocks 



(1) 



where we defined a two component field vector £ = 
(f-l-,f_)^, a two-component vector of spin polariza- 
tions S = {^^ags, \/nag}iY' , and a two-component vec- 
tor of optical polarizations P = {\pnage.j^,\fnageJ)^ ■ 
Here K_ = az-§^ and G_ = 9\/n l2x2, where g = 
pyijj/2h£^ is the coupling constant of £± to the tran- 
sition \g) — |e) with dipole matrix element p and n is the 
atomic number density. Furthermore 5_ = diagjcJ,,, 5^}, 

o — (^s+ ^h+ 

- ~ ^- nu-^ 

Finally Fp are Langevin noise forces neccessary for pre- 
serving the commutation relations of the decaying vari- 
ables P. For an exponential decay these operators are 
(5-correlated in time, i.e. {FA{t)FB(t')) = DABS{t — t'), 
where the diffusion coefficients Dab are proportional to 
the population of the excited states and can be calcu- 
lated by means of the dissipation-fiuctuation theorem 
[26] . With the help of a generalized Morris-Shore trans- 
formation [27], one can show that in linear response this 
system has two adiabatic eigensolution that are decou- 
pled from the excited states. If we choose the control 



, r = diag{7ge+ - iA+,7ge_ - iA„}. 



fields such that £2 = -^^(1 + ia^), where ax is a Pauh 
matrix, the eigensolution of the tripod systems, the dark- 
state polaritons (DSP), are 



^+ — cos9£+ j= sinfl {a-gs — i<7gh) , (2) 

^- = cos9£^ H — -= sind {i<Tgs — (Tgh) , (3) 

v2 

where tan^ 9 = g^n/n^ is the mixing angle between light 
and matter excitation. We emphasize that in contrast 
to ordinary stationary light where only one indepen- 
dent dark polariton exists [27] we here have two inde- 
pendent dark eigensolutions. We now introduce a small 
two-photon detuning of opposite sign 5 = Sg ~ —Sh 
for the transitions \g) — \s} and \g} — |/i), respectively, 
with |(5| <C ^o/t- Assuming one-photon resonance, i.e. 
A± = 0, setting 7ge± = 7, and transforming cq. (1) 
to the polariton basis and adiabatically eliminating all 



variables except for the DSPs, we arrive at 
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iiabsfg sin ^ ( Cr^^ CTy ) * + -F*- (4) 



Eq. (4) represents a Dirac equation for a spinor ^ 
{^+, 'J-) with effective "speed of hght" and mass 
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The imaginary quadratic term in (4) results from non- 
adiabatic couplings to decaying states, F* is the corre- 
sponding Langevin-noise operator, iabs ~ cj/g'^n is the 
resonant absorption length in absence of EIT, az and 
ay denote Pauli matrices. From eq. (4) we expect a 
Dirac cone structure of the dispersion relation as seen 
in Fig. 2. The interaction of the control fields with the 
atoms results in two Autlcr-Townes states represented 
by the horizontal black dashed lines. The dispersion of 
the free probe fields is given by the almost vertical black 
dashed lines. Switching on the interaction, the red dotted 
dispersion branches are formed. These correspond to the 
two components of the slow-light spinor and are in a very 
good approximation described by the expected Dirac dis- 
persion with c* = Vg around k ~ (orange dashed line). 
Finally inducing a small two-photon detuning S results 
in a finite mass m*, as can be seen by the finite splitting 
of the blue (solid) line. 

To confirm the effective model, we numerically inte- 
grated the full set of MBE for the scheme of Fig. Id in 
time for the case of m* = 0. The results showed very 
good agreement with the analytic prediction. 

We now exploit the property of a locally adjustable 
mass of the SSP. If at a certain point along the z-axis 
the mass changes its sign, i.e. ttiq — >■ —'m-Oi one can eas- 
ily show, that there exists a localized mid-gap (i.e. zero 
energy) solution around the mass jump, as by squaring 
the Dirac Hamiltonian leads to a pair of supersymmctric 
Hamiltonians corresponding to two independent parti- 
cles, described by a Schrodinger equations with a Dirac-J 
potential. The zero-energy bound state is exponentially 
localized with a localization length proportional to the 
effective Compton length X^ = h/{m*c*). This is indi- 
cated in Fig. 3. 

Thus one can ask how the state will look like when 
the number of mass jumps and their size is increased? 
This leads to the ID random-mass Dirac model: Given a 
spatially random mass m{z)^ with (ft = c = 1) 



m{z)m{y) = 2V5{z ^ y), m{z) = 0, 



(6) 



where F parameterizes the strength of the disorder, it 
was shown [4, 5] that the density of states at zero energy 
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FIG. 3. (Color online) Zero energy bound state at a mass 
jump 771(2) = msgn{z). a = Xc — h/{m*c') is the effective 
Compton length. 
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FIG. 4. (Color online) Correlation function \I{z)I(0)\ with 
disorder width a = 0.0 17 averaged over 50 realizations, or 
large times, the numerical result show very good agreement 
with the theoretical expected power-law (black dashed line), 
with the only fit-parameter being the ordinate intersection. 
The inlet shows the intensity I(t) spatially integrated over a 
region ±3Lp. The curve is normalized to the initial intensity 
I{to, z = 0). The fitted curve (blue crosses) corresponds to the 
theoretically predicted effective decay rate with L* ~ \/2I/p. 



is diverging, corresponding to a localized state. The cor- 
responding wave function 4'(z) exposes however unusual 
density correlations. Contrary to expectations there is 
not an exponential localization as in the case of Fig. 3 
but one finds for Fz ^ 1: 



|/(z)J(0)| 



3/2 



(7) 



where I{z) = |*(z)p. A spatially varying mass (de- 
tuning S) can be achieved by using fluctuating magnetic 
fields or by applying speckle patterns as done for the 
demonstration of Anderson localization [28] . 

In experiments one has a finite correlation length ^ of 
the disorder which we use as a discretization length for 
the simulations. We assume that the two-photon detun- 
ing has a Gaussian distribution with 6 = and 5'^ = a^ 
for all points in space. 

In Fig. 4 numerical results of the density correlation 



(7) are shown obtained from the full MBE. For large 
times the correlation (red solid line) approaches the pre- 
dicted power-law behavior, shown by the black dashed 
line. 

The finite two-photon detuning also leads to some 
small losses, seen in the inlet of Fig. 4. Two regimes are 
apparent: The first rapid loss is due to imperfect match- 
ing of the initial wave function to the zero-energy eigcn- 
statc. The smaller losses for large times can be attributed 
to non-adiabatic couplings of the SSP to decaying states, 
resulting in an effective loss rate jcS = I'^iabsi'g/i*, 
L* being some characteristic length scale of the localized 
state. These losses are negligible as long as jctfT <C 1, 
where T = L/vg is the typical time scale of the ex- 
periment. Furthermore (7) requires TL ^ 1, in or- 
der to observe the power-law decay. This translates to 
L » Wgnkinks/c^ sin'' 9 = Leon-, wlicrc TT-kinks is the num- 
ber density of mass jumps. Both equations finally lead to 

L/Labs = OD » max{gv^]Vk~^,(7rL/i,)'} • Here 
-^kinks = 'T'kinksi is the number of kinks in the length L 
and OD denotes the optical depth of the medium. For 
white noise the mean free path Amcan = "-kinks '^^^ ^"^ 
approximated by the disorder correlation length ^, i.e. 
Amoan ~ C- Taking ^ « O.S/im, L « 4 mm [28], and 
estimating A'kinks = i/^f ~ iC yields a required optical 
depth of OD > 100. This high optical depth is however 
achievable in systems such as atoms trapped in a hollow 
core photonic bandgap fiber [29] or cold atoms trapped 
in the evanescent fields of an ultra-thin optical fiber [30] . 

An estimate of the minimal build-up time of the cor- 
relations, defined as Tcorr = Lcoi-r/vg yields Tcorr ~ 
2.5 - 25fis. Here we used a w O.I7 == (27r) 0.6MHz (for 
^^Rb), Wg « 10 — lOOm/s [3] and the correlation length 
^ « 0.3/im from above. This time is far shorter than 
the achievable lifetime of the dark-state polaritons, typi- 
cally on the order of 100/is [6] to several ms [11, 12], thus 
enabling experimental observation. 

A possible experimental protocol to observe the pro- 
posed effect is to store an initial wave packet in the 
atomic coherences with existing storing techniques for 
slow-light polaritons [6] or by coherent RF-transitions 
[31]. Subsequently one applies the disorder, e.g. by a 
laser-speckle induced ac-Stark shifts and reads out the 
polaritons by applying all four control fields. The polari- 
ton states then evolve according to the above presented 
model. After a certain time the disorder is removed and 
the final state is read out by means of usual slow- light. 
Looking at the time profile of the intensity of the probe 
fields arriving at the detector one can reconstruct the 
spatial intensity profile of the localized state. 

In summary we have shown that the coupling of two 
counter-propagating light fields to an optical thick en- 
semble of atoms driven by multiple drive fields in a 
tripod-like linkage pattern leads to the formation of 
spinor-like eigcnsolutions which are immune to sponta- 



neous emission losses. These spinor dark-state polaritons 
obey a ID Dirac equation with effective speed c* given by 
the group velocity of slow- light in FIT media and with an 
effective mass m* determined by a small two-photon de- 
tuning. This allows to experimentally study a number of 
interesting phenomena of relativistic quantum dynamics. 
For example the possibility to tune the effective mass can 
be used to observe the unusual localization transition pre- 
dicted for the random-mass Dirac model, which appears 
in the context of several interesting condensed matter 
systems with disorder. Furthermore adding interactions, 
e.g. by means of resonantly enhanced Kerr nonlinearities 
[32] , important relativistic many-body models such as the 
bosonic counterpart of the Thirring model [33] should be 
experimentally accessible. 
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